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ABSTRACT 



We present new high-resolution observations of Sagittarius A* at wavelengths 
of 17.4 to 23.8 cm with the Very Large Array in A configuration with the Pie 
Town Very Long Baseline Array antenna. We use the measured sizes to calibrate 
the interstellar scattering law and find that the major axis size of the scattering 
law is smaller by ~ 6% than previous estimates. Using the new scattering law, 
we are able to determine the intrinsic size of Sgr A* at wavelengths from 0.35 
cm to 6 cm using existing results from the VLBA. The new law increases the 
intrinsic size at 0.7 cm by ~ 20% and < 5% at 0.35 cm. The intrinsic size is 13^3 
Schwarzschild radii at 0.35 cm and is proportional to A 7 , where 7 is in the range 
1.3 to 1.7. 

Subject headings: Galaxy: center — galaxies: active — techniques: interferomet- 
ric — scattering 



1. Introduction 

Imaging the radio emitting region surrounding the massive black hole in the Galactic 
Center, Sagittarius A*, has been a goal since its discovery (Balick & Brown 1974). Turbulent 
electrons along the line of sight to Sgr A*, however, scatter radio wavelength photons and 
produce an image that is an elliptical Gaussian with a major axis size of ~ 0.5 arcsec at 20 cm 
and a A 2 dependence (Backer 1978). Separating the effects of the small intrinsic source from 
the effects of scattering has required observations at short wavelengths, careful calibration, 
and the use of closure amplitude techniques, which reduce sensitivity but remove uncertainty 
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due to calibration error (e.g., Bower et al. 2004; Shen et al. 2005). These efforts have recently 
resulted in the first robust detections of intrinsic structure in Sgr A* at wavelengths of 1.3 
cm, 0.7 cm and 0.35 cm. The intrinsic source has a size that scales with A 11 or A 16 to a 
minimum of ~ 10 Schwarzschild radii at 0.35 cm (assuming = 4 x 1O 6 M and d = 8 
kpc ; Ghez et al. (2005); Eisenhauer et al. (2003)). 

These detections of the intrinsic size of Sgr A* have a number of consequences. The 
brightness temperature of 10 10 K strongly excludes advection dominated accretion flows 
(ADAFs; Narayan et al. 1998) and Bondi-Hoyle accretion models (Melia 1994). These size 
measurements, however, cannot differentiate between jet models (Falcke et al. 1993), generic 
radiatively inefficient accretion flows (Quataert & Gruzinov 2000), and hybrids of these 
models (Yuan et al. 2002). This limitation is principally due to the limited sensitivity in the 
minor axis size of the scattering ellipse. Coupled with measurements of the proper motion 
of Sgr A* (Reid & Brunthaler 2004), the assumption that the black hole is smaller than the 
emission region implies a lower limit to the mass density of the black hole ~ 10 5 M Q AU~ 3 , 
which strongly excludes alternative models for dark mass objects. 

The scattering medium itself is a system of intense interest (Lazio & Cordes 1998; 
Bower et al. 2001; Goldreich & Sridhar 2006). The A 2 size dependence of Sgr A* is a strict 
consequence of the strong scattering and the short projected baselines at the distances of 
the scattering medium (Narayan & Goodman 1989). 

We present here detailed measurements of the scattering properties at wavelengths that 
range from 17.4 to 23.8 cm using the Very Large Array and the Pie Town Very Long Baseline 
Array antenna (§2). The addition of the PT antenna to the VLA A configuration improves 
the East-West resolution by a factor of two. The resulting scattering law is smaller by 
~ 6% than previous estimates. The independent estimate of the scattering law enables us to 
measure the intrinsic size of Sgr A* at a wavelength as long as 6 cm (§3). We discuss these 
results in (§4). 

We wish to clarify here the numerous axes and steps involved in translating observations 
of Sgr A* into a measurement of the intrinsic size. VLA observations of Sgr A* are obtained 
with a synthesized beam that is extended North-South. Deconvolution of the observed 
image with the synthesized beam gives the apparent, scatter- broadened image of Sgr A*. 
This image is predominantly a two-dimensional Gaussian with the major axis oriented ~ 80 
degrees East of North. Throughout this paper when we refer to the major and minor axes, 
we refer to the orientation of the scattering Gaussian. Finally, to obtain the intrinsic image, 
we deconvolve the apparent image with a model of the scatter-broadened image, which is 
determined from long wavelength apparent sizes. 



- 3- 



2. Observations and Analysis 

We obtained new observations of Sgr A* with the Very Large Array plus the Pie Town 
VLBA antenna. The VLA was in the A configuration. Observations were made on 1 and 4 
October 2004 in eight separate bands centered at 25.2, 23.8, 23.2, 21.9, 20.9, 19.8, 18.0, and 
17.5 cm. Each band had 12.5 MHz bandwidth with 15 channels. Results at 25.2 and 21.9 
cm were corrupted by interference and we do not consider these data any further. 

We calibrated the absolute flux density with observations of 3C 286. Corrections for 
atmospheric and instrumental amplitude and phase fluctuations were made through self- 
calibration of frequent observations of the compact source Jl 744-312. We imaged Sgr A* 
using baselines longer than 50k\ and uniform weighting to suppress large-scale structure 
(Figure 1). 

The presence of a radio transient with a flux density of ~ 30 mJy and a resolved 
morphology at a separation of 2.7" South of Sgr A* precluded modeling in the visibility and 
closure amplitude domains (Bower et al. 2005). We previously showed that fitting in the 
image domain provides results that are equivalent to fitting in the closure amplitude domain 
at centimeter wavelengths, in the case where the difficulty of calibration and poor telescope 
performance are less critical (Bower et al. 2001, 2004). Long wavelength data obtained from 
the VLA meet these criteria better than any other data. Accordingly, we fit sizes to Sgr A* 
in the image plane with a region that excluded most of the transient flux density. The effect 
of the transient is primarily on the accuracy of the size in the minor axis of the scattering 
angle. The synthesized beam ranges from 1.69 x 0.56 arcsec at 17.5 cm to 2.36 x 0.98 arcsec 
at 23.8 cm, oriented roughly North-South. Measured sizes ranged from 1.71 x 0.67 arcsec at 
17.5 cm and 2.43 x 1.22 arcsec at 23.8 cm. Sgr A* is clearly resolved in both axes but with 
considerably more accuracy in the East- West axis than an in the North-South axis. Fitting a 
point source to the data produced very poor quality fits, while fitting an elliptical Gaussian 
produced a residual image with no obvious systematic errors and an rms ~ 2.5 mJy/beam 
(Figure 1). This rms is a few times the rms ~ 0.9 mJy/beam determined far away from Sgr 
A*, possibly due to the presence of confusing emission around Sgr A*. 

We deconvolved the measured Gaussian sizes with the synthesized beam sizes to deter- 
mine the true source parameters: total flux density, major axis, minor axis, and position 
angle. We determined errors in the parameters by calculating \ 2 f° r a grid in the parame- 
ters surrounding the best-fit value. These errors are the formal uncertainy in the parameters 
and do not reflect the systematic errors, which we discuss below. Results are tabulated in 
Table 1. 

In the case of the minor axis, there is a clear trend of decreasing size with decreasing 
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wavelength, which indicates the presence of systematic errors. Marginal resolution of Sgr A* 
in the North-South axis and confusion from the presence of the radio transient due South of 
Sgr A* are the likely causes of this effect. The major axis size, however, is not affected by 
the transient and only weakly distorted by changes in the position angle; the best-fit solution 
for major axis size changes by only 1% with a 10° change in position angle. We conclude 
that the results determined from the VLA+PT result are accurate in the major axis but not 
in the minor axis. 

We experimented with a range of imaging parameters to explore systematic effects on 
the deconvolved size of Sgr A*. Weighting with a robustness parameter of 0, using super- 
uniform weighting, and changing the minimum (u, v) distance from 40 to 100 kX changed the 
deconvolved major axis size by no more than 3%. Since the results are strongly dependent 
on the PT antenna, we dropped random groups of 5 baselines associated with PT, producing 
1% changes in the deconvolved size. These errors are comparable to those found for other 
sources through VLA observations (e.g., Trotter et al. 1998). Thus, our results are influenced 
by systematic imaging effects at a level of a factor of no more than a few. As we discuss 
below, a factor of a few is consistent with the scatter in the measurements. 

3. The Scattering Size and The Intrinsic Size of Sgr A* 

In Figure 2, we plot the measured size of Sgr A* from the VLA+PT observations and 
from Very Long Baseline Array results from Bower et al. (2004) at wavelengths of 6.01 cm 
to 0.67 cm. We also include the VLBA result from Shen et al. (2005) at 0.35 cm. The sizes 
are plotted normalized by wavelength squared. 

We fit a power law of the form A 2 to the major axis sizes using the new VLA+PT results 
at 17.4 to 23.8 cm wavelength. The best-fit value to the normalization of the scattering law 
is 1.309 ± 0.017 mas/cm 2 . The errors in these values are determined from the scatter in the 
measurements. The best-fit scattering values are plotted as a straight line in Figure 2. The 
major axis normalization is ~ 6% smaller than previous estimates. None of our measurements 
are consistent with the previously determined major axis normalization of 1.39 mas/cm 2 . 
The best estimates of the minor axis scattering size and position angle remain the results 
determined previously from VLBA observations at wavelengths between 2 cm and 6 cm 
(0.64+°;^ mas/cm 2 and 78+°;^ deg; Bower et al. 2004). 

The scatter in the VLA+PT major axis sizes is much larger than the expectation of the 
statistical errors for individual points. The reduced xl ~ 7 for the major axis, indicating that 
there are additional sources of error in the measurement of the size that we are not including. 
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Dropping either the two highest frequency or two lowest frequency VLA+PT sizes did not 
significantly affect the reduced x 2 or t ne best-fit scattering law. We also explored the effect 
of the inclusion of the shorter wavelength VLBA results on the major axis scattering size. 
Inclusion of the 6 cm size decreases the reduced xl slightly but does not affect the scattering 
size significantly. Inclusion of VLBA results at wavelengths of 3.6 cm and shorter, however, 
leads to a significant increase in xl to 18- That is, the major axis sizes at A <= 3.6 cm are 
not consistent with the longer wavelength sizes and a A 2 scattering law. This inconsistency 
holds if we scale the VLA+PT errors by a factor as large as 5, which reveals xl — 1-8- 

The A 2 dependence of the scattering law is strongly favored for theoretical reasons. The 
maximum baseline length projected to the scattering region is b pro j = D scatterin g/ Ds gr A x 
bmax ~ lkm, where D Sgr A = 8 kpc is the distance to Sgr A*, D scattering ps 100 pc is 
the distance of the scattering region from Sgr A*, and b max ps 70 km is the maximum 
baseline between the VLA and PT. b pro j is substantially smaller than the expected and 
measured inner scales (bi nner ) for the power spectrum of turbulent fluctuations (10 2 to 10 5 ' 5 
km; Wilkinson et al. (1994)). This result holds for the long VLBA baselines involved in 
imaging at shorter wavelengths, as well, where b pro j ~ 25 km at 0.7 cm wavelength. For the 
case of bp ro j << b inner , the resulting image is very heavily averaged and must be Gaussian 
in shape with size oc A 2 (Narayan & Goodman 1989). 

This expectation of strong scattering is supported by previous measurements of the 
shape of the image of Sgr A*. Bower et al. (2004) showed that fitting the closure amplitudes 
of Sgr A* at 0.7 cm with a functional form for the visibilities of oc e~ b</J 2) , where b is the 
baseline length, revealed (5 = 4.00 ±0.03. That is, the best evidence indicates that the image 
of Sgr A* is a Gaussian. Following scattering theory, where the size is proportional to \ a , 
then a = /3-2 = 2.00 ± 0.03 (Narayan & Goodman 1989). For the case of the VLA+PT 
data, we find a = 1.98 ±0.11, which is consistent with the expectation of a — 2. 

Nevertheless, if we fit a single power-law to all of the VLA+PT and VLBA data, we find 
a = 1.96 ± 0.01 with a xl = 5-5. So, without the constraint of a A 2 size for the scattering 
law at some wavelength, the evidence for resolving an intrinsic size becomes marginal at 
any wavelength. A fit of the size proportional to Va 2 A 4 + 6 2 A 27 produces an apparent size 
very similar to that of the single index power-law fit but is not sufficiently constrained to 
set reasonable limits on the parameters. If we fix 7 and search for a and b, we find that 
X 2 for 7 = 1 is four times the value for 7 unconstrained, indicating that 7 = 1 is strongly 
excluded. Without the assumption that the second term is negligible for wavelengths longer 
than 6 cm, therefore, we cannot determine the scattering law or the intrinsic size of Sgr A*. 

A final caveat is required. The scattering medium is dynamic. The minimum time scale 
for a change in the medium is the refractive time scale, which is 0. by {v /100km/ s){\/ lcm) 2 
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for Sgr A*, where v is the velocity of the scattering material relative to Sgr A* and the Sun 
(Narayan & Goodman 1989). The data presented here were obtained in a span of roughly a 
decade. The long-wavelength scattering properties are very unlikely to change on this time 
scale. However, at wavelengths as long as 4 cm, the refractive time scale is < 10 year. The 
many observations at 0.7 cm in this period, however, appear to produce a source of stable 
size, despite a refractive time scale less than one year. We conclude it is unlikely that the 
scattering size has changed significantly over this period. 

With the assumption that the scattering law is determined accurately at wavelengths 
longer than 17 cm, we can determine the intrinsic size. We subtract in quadrature the 
scattering law size from the measured size (Table 2, Figure 3). We compute the results for 
the best-fit major axis scattering law {b sc = 1.31 mas/cm 2 ), and ±3a of the best-fit value. 
For the best-fit case, the intrinsic size is accurately determined from 0.35 cm to 3.6 cm. Over 
this range, the intrinsic size is well-fit by a power-law A 7 , where 7 = 1.6 ±0.1. For the smaller 
scattering size, we find a steeper power-law and measure the intrinsic size from 0.35 cm to 6 
cm. For the larger scattering size, we cannot measure the intrinsic size at wavelengths longer 
than 1.3 cm and find a shallower power-law index of 7 = 1.3 ± 0.2. 

If the intrinsic size power-law extends to A ~ 20 cm, the contribution of the intrinsic 
size results in an increase of the measured angular sizes by ~ 1.5%. This is comparable to 
the error in the major axis scattering law and, therefore, negligible. 

4. Discussion 

We have measured the intrinsic size of Sgr A* from 0.35 to 3.6 cm. At short wavelengths, 
the result is consistent with the conclusions of recent efforts by Bower et al. (2004) and Shen 
et al. (2005). The size of Sgr A* at 0.35 cm is 13.3l|;Ji2 a , where R s = 1.2 x 10 12 cm is the 
Schwarzschild radius for = 4 x 1O 6 M and the Galactic Center distance d — 8 kpc. 
This compact size confirms tight restrictions on accretion models and black hole alternatives 
previously claimed and stated in §1. 

The wavelength dependence of the source size agrees with that found by Bower et al. 
(2004) and is steeper than that found by Shen et al. (2005), who found 7 ~ 1.1. The steeper 
dependence indicates that the brightness temperature decreases as A -1 , assuming that the 
size in the second dimension is proportional to the major axis size. The peak brightness 
temperature at 0.35 cm is ~ 10 10 K for a flux density of 1 Jy. The power-law dependence of 
the size as a function of wavelength indicates a stratified, smoothly varying emission region. 

Detailed jet models for Sgr A* predict 7 ~ 1 (Falcke & Markoff 2000). Generalized 
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jet models, however, allow a range of 7, depending on the details of the magnetic field and 
particle energy density distributions (e.g., Konigl 1981). A jet with B oc r -1 , electron density 
decreasing as r" 1 , and optically thin power-law index of 1 will show a size oc A 1 ' 4 . 

Yuan et al. (2006) model the size of Sgr A* for a radiatively inefficient accretion flow. 
They fit sizes at 0.35 cm and 0.7 cm that are fit with a power law of index 7 = 1.1. Variations 
in the nonthermal electron distribution or deviations from equipartition, however, could alter 
7 in their model. 

The critical remaining observational goals for understanding the image of the radio emis- 
sion of Sgr A* are a measurement of the two-dimensional size and detection of structural 
variability. The simple one-dimensional deconvolution that we have performed here only 
gives schematic information on the size of Sgr A*. With a more accurate two-dimensional 
scattering model, future analysis will directly compare the observed image with non-Gaussian 
source models convolved with the imaging and scattering constraints. Astrometric observa- 
tions may indicate a shift in the centroid of the image with frequency. A heterogeneous jet 
will exhibit such a shift due to changing location of the optically thick surface of the source 
(Konigl 1981). 

At mm and submm wavelengths, the gravity of the black hole will distort the image 
(Falcke et al. 2000; Broderick & Loeb 2005). Detailed knowledge of the shape of the longer 
wavelength image will permit a more accurate characterization of light-bending effects in the 
actual image. Ultimately, these images will provide one of the strongest tests of the existence 
and characterization of black holes. 

The National Radio Astronomy Observatory is a facility of the National Science Foun- 
dation operated under cooperative agreement by Associated Universities, Inc. 
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Fig. 1. — (left) Image of Sgr A* at 17.5 cm from the VLA+PT observations. Sgr A* is 
the bright source at the center of the image. The source to the South of Sgr A* is the 
transient. The synthesized beam size is indicated in the lower part of the image. Contours 
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begin at 5 mJy/beam and increase by factor of y(2) to half the peak intensity of Sgr 
A*. The synthesized beam size is 1.69 x 0.56 arcsec. (center) Residual image of Sgr A* 
after subtraction of best-fit point source. This shows that the source is clearly resolved 
East- West and at best marginally resolved North-South, (right) Residual image of Sgr A* 
after subtraction of best-fit Gaussian source. The source is well-modeled to an rms of 2.5 
mJy/beam. 
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Fig. 2. — Measured major axis size as a function of wavelength. Triangles are VLBA mea- 
surements determined through closure amplitude analysis from Bower et al. (2004) and Shen 
et al. (2005). Squares are the new VLA+PT measurements. The major axis sizes have been 
normalized by A 2 . The solid red line is the best-fit scattering value determined from the 
VLA+PT data alone. Dotted red lines are ±3cr of the best-fit scattering law. 
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Fig. 3. — Intrinsic size in Schwarzschild radii for the East- West dimension using three dif- 
ferent estimates of the major axis scattering law. We assume a 4 x 1O 6 M black hole at a 
distance of 8 kpc. In the central panel, we show results for the best-fit scattering law. In 
the left and right panels, we show the results for the — 3cr and +3a, respectively. The solid 
red lines are the best fit curves for size oc A 7 . 
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Table 1. Apparent and Deconvolved Sizes of Sagittarius A* from VLA+PT Observations 

of October 2004 



Apparent Size Deconvolved Size 

A Synth. Beam b maj b min PA b maj b min PA S 

(cm) (mas X mas,deg) (mas) (mas) (deg) (mas) (mas) (dcg) (mjy) 



23.8 
23.2 
20.9 
19.8 
18.0 
17.5 



2357 x 979, 0.1 
2191 x 714, 4.8 
1996 x 646, 7.6 
1840 x 581, 10.4 
1737 x 564, 7.2 
1689 x 556, 6.8 



2433.5to'o 1222.6±H -0.21 



2049.4Z 



2252.1+j5 
f6 
-6 

1877.8^5 
1758.612 
1706.4+j 



3.1 

985.5lg° 

849.2I2I 
,+2:2 
-2.2 

695.9lJo 
-1.0 
-1.0 



751. r 



667.6 



1+0.14 
-0.14 
- 00+O.IO 

°-° a -o.w 

q n 6 +0 07 
y ' uo -0.07 
11.92+ 06 



-0.06 
1-0.06 
-0.06 
7 71 +0.06 



8.23 



734.7tg;5 
697.lt*;* 
601.8+Jg 

i+6'.3 
'-8.0 

430.8l 2 -g 

7+2:4 

-3.0 



522.8 



391.27 



602A+f g ° 9 
496.3tS;j 

+27.2 
19.1 
+31.2 
-21.1 

236.7±ilt 
+21.6 

15.5 



397.3 
306.6 



206.3 



97.7 
76.5 
66.1 
70.5 
74.9 
74.5 



+ 1.1 

7.5 
+2.7 
2.1 
3.9 
2.9 
+3.7 
2.7 
+ 1.8 
1.3 
+2.1 
1.5 



471. 
490. 
432. 
441. 
505. 
511. 



1.6 

1.6 
1.5 

-1.9 

+ 1.3 
1.3 
1.3 
1.3 

+0.7 
0.7 

+0.7 
0.7 



Note. — The quoted errors in the peak flux density S do not include the ~ 10% error in absolute flux density. 



Table 2. Intrinsic Size of Sagittarius A* 



b sc = 1.26 mas/cm 2 b sc = 1.31 mas/cm 2 b sc = 1.36 mas/cm 2 
A Size Size Size 

(cm) (mas) (mas) (mas) 



0.35 


f) 142+0-062 
U.1^Z_ 025 


1 Qfi+0.069 
U.lOU_ 


n 1 OC+0.076 
U.lZO_o.040 


0.67 


n qoc;+o.i57 


362 +ai9 ° 

u - ouz -0.080 


0.3241°;^ 


0.69 


n Q7Q+o.iio 
u.o i o_ 077 


n oqi +0.152 

u.oo±_ 120 


0.2801°;^ 


0.76 


n 5n5+ - 177 

u.ouo_ 128 


n 4fii + - 221 

U.4Ol_ 172 


0.410;°;I 2 


1.35 


i 9QC+0.369 
i.zyo_ 174 




8Q7 +0 ' 767 

u - oy ' -0.571 


1.95 


9 fi98+ - 408 
Z.UZO_ 304 


1.50018;^ 


u.duo_ 562 


3.56 


7 ZLQ1+ 1 - 524 

' - 4yi -0.281 


5.9531?:$ 


o 7 qp;+5.280 
o. i oo g 732 


6.01 


16 427+ 2 - 981 

iD - 4Z ' -16.410 


10.094l? ^3 


o ooo+ 19 - 408 

u.uuu_ 000 



Note. — The intrinsic size is computed for the best-fit value to the major axis scattering 
size {b sc = 1.31 mas/cm 2 ) and for ±3<r values. 



